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On  The  Unsteady-Motion  Theory  of  Magnetic  Forces  for  Magiev 

by 

S.  S.  Chen,  S.  Zhu,  and  Y.  Cai 


Abstract 


Motion-dependent  magnetic  forces  are  the  key  elements  in  the  study  of 
magnetically  levitated  vehicle  (magiev)  system  dynamics.  In  the  past,  most  maglev- 
system  designs  were  based  on  a  quasisteady-motion  theory  of  magnetic  forces.  This 
report  presents  an  experimental  and  analytical  study  that  will  enhance  our 
understanding  of  the  role  of  unsteady-motion-dependent  magnetic  forces  and 
demonstrate  an  experimental  technique  that  can  be  used  to  measure  those  unsteady 
magnetic  forces  directly.  The  experimental  technique  provides  a  useful  tool  to 
measure  motion-dependent  magnetic  forces  for  the  prediction  and  control  of  magiev 
systems. 


1  Introduction 


The  dynamic  response  of  magnetically  levitated  vehicle  (magiev)  systems  is 
very  important  for  the  safety,  ride  quality,  guideway  design,  and  cost  of  these 
systems.  The  key  elements  in  the  study  of  magiev  system  dynamics  are  magnetic 
forces.  In  the  past,  many  studies  have  been  performed  to  help  us  understand  and 
predict  magnetic  forces  for  differing  cases.  It  appears  that  most  maglev-system 
designs  are  based  on  the  magnetic  forces  calculated  or  measured  with  the 
quasisteady-motion  theory. 

Some  limited  studies  of  magnetic  forces  that  consider  the  effects  of  unsteady 
motion  have  been  performed.  Magnetic  damping  was  considered  by  Fujiwara 
(1980),  Iwamoto  et  al.  (1974),  Moon  (1977),  and  Yamada  et  al.  (1974).  In  these 
investigations,  the  response  of  a  system  supported  by  or  subjected  to  magnetic 
forces  was  measured  or  analyzed.  From  the  response  of  the  system,  the  magnetic 
damping  was  calculated.  This  can  be  considered  an  indirect  method  to  quantify 
magnetic  forces.  Another  method  is  to  calculate  the  magnetic  damping  by  using 
eddy  current  losses  (Saitoh  et  al.  1992)  or  the  magnetic  forces  due  to  arbitrary 
motion  (Davis  and  Wilkie  1971).  These  are  approximate  solutions  that  use  a  d:rect 
method  that  is  based  on  quasisteady  motion. 

The  indirect  method  has  certain  limitations.  Magnetic  coupling,  including 
magnetic  damping  and  stiffness,  is  difficult  to  quantify.  A  direct  method  will  be 
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more  useful  for  identifying  each  element  of  the  magnetic  forces.  In  addition,  it  is 
not  clear  whether  the  quasisteady-motion  theory  is  adequate  for  the  development  of 
a  stability  theory  for  maglev  systems.  In  some  cases,  the  unsteady  magnetic  forces 
are  important  and  the  quasisteady-motion  theory  appears  to  be  inadequate.  In 
developing  future  maglev  systems,  at  least  some  scoping  calculations  on  the  effects 
of  unsteady  magnetic  forces  must  be  performed  to  avoid  unexpected  failure  or 
excessive  operating  and  maintenance  costs. 

This  report  presents  a  general  approach  to  the  calculation  of  motion-dependent 
magnetic  forces  and  an  experiment,  based  on  the  unsteady-motion  theory,  to 
measure  magnetic  forces.  The  objectives  of  this  study  are  to 

•  Understand  the  role  of  motion-dependent  magnetic  forces. 

•  Provide  some  theoretical  considerations  of  analytical/numerical  methods 
for  calculating  unsteady  magnetic  forces. 

•  Develop  and  demonstrate  an  experimental  technique  that  can  be  used  to 
measure  motion-dependent  magnetic  forces. 

This  study,  for  the  first  time,  presents  a  direct  method,  based  on  the  unsteady- 
motion  theory,  to  measure  motion-dependent  magnetic  forces. 


2  Motion-Dependent  Magnetic  Forces _ 

Magnetic  forces  must  be  considered  in  any  analysis  of  vehicle  dynamics, 
guideway  structural  and  fastening  design,  and  prediction  of  ride  comfort  (Chen  et 
al.  1992).  These  forces  are  considered  on  the  basis  of  how  they  relate  to  vehicle 
stability.  In  a  vehicle,  many  magnets  or  coils  are  used  in  the  levitation,  guidance, 
and  propulsion  systems.  Without  loss  of  generality,  a  single  magnet  or  a  single  coil 
is  used  here  as  an  example  to  demonstrate  the  role  of  motion-dependent  magnetic 
forces. 

Consider  a  coil  or  a  magnet  with  three  degrees  of  freedom,  three  translations  in 
the  propulsion  (x),  guidance  (y),  and  lift  (z)  directions  as  shown  in  Fig.  1.  Let  U  be 
the  vector  consisting  of  the  three  motion  components,  i.e., 


U1 

U2 

•  =  ■ 

uy 

u3 

uz 

(1) 
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Velocity  and  acceleration  are  given  by 


and 


(2) 


a2u 

at2' 


The  magnetic  forces  in  the  three  directions  can  be  decomposed  into  two  parts, 
steady  magnetic  forces  and  motion-dependent  magnetic  forces. 


2.1  Steady  Magnetic  Forces 

When  the  coil  or  magnet  is  stationary,  the  three  components  of  the  steady 
magnetic  forces  fxo>  fy0,  and  f20  can  be  written 


[flo] 

fxo 

II 

«2o 

II 

fyo 

f3o. 

fzo 

The  steady  magnetic  forces  are  independent  of  the  motion  of  the  magnet  or  coil. 


2.2  Motion-Dependent  Magnetic  Forces 


When  the  coil  or  magnet  is  moving,  the  motion-dependent  magnetic  forces  fx, 
fy,  and  fz  can  be  written 


fl 

fx 

F  =  - 

h 

►  =  * 

fy 

f3 

fz 

(4) 


The  explicit  form  of  the  motion-dependent  magnetic  forces  is  not  available.  It  can 
be  written  as 

3 

fi  =  Z(miiaj+cij'ij+kijuj)’  i=1>2>3> 

j=l 


(5) 
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where  mjj,  cjj,  and  kjj  are  magnetic  mass,  damping,  and  stiffness  coefficients.  Those 
coefficients  are  functions  of  the  system  parameters  and,  in  general,  they  are  also 
functions  of  uj,  uj,  and  Uj. 

The  motion-dependent  magnetic  forces  have  been  considered  from  the 
viewpoint  of  two  theories. 


2.3  Quasisteady-Motion  Theory 

The  magnetic  forces  acting  on  an  oscillating  magnet  are  equal,  at  any  instant 
in  time,  to  those  of  the  same  magnet  moving  with  a  constant  velocity  and  with 
specific  clearances  equal  to  the  actual  instantaneous  values.  The  magnetic  forces 
depend  on  the  deviation  from  a  reference  state  of  speed  and  clearance,  i.e.,  the 
motion-dependent  magnetic  forces  depend  only  on  U,  but  not  on  U  and  U.  In  this 
case,  mjj  and  cij  are  zero.  Only  magnetic  stiffnesses  are  considered. 

2.4  Unsteady-Motion  Theory 

„  The  magnetic  forces  acting  on  an  oscillating  magnet  will  depend  on  U,  U,  and 
U,  as  shown  in  Eq.  5.  The  magnetic  force,  based  on  the  unsteady-motion  theory, 
can  be  obtained  by  calculating  or  measuring  the  magnetic  force  acting  on  the 
magnet  oscillating  in  the  magnetic  field. 


3  Analysis  and  Experimental  Methods _ 

Two  methods  can  be  used  to  study  the  motion-dependent  magnetic  forces,  an 
indirect  method  and  a  direct  method. 


3.1  Indirect  Method 

A  magnet  or  coil  that  is  levitated  magnetically  can  be  analyzed  theoretically  or 
numerically  to  predict  its  response.  From  the  time  histories  of  the  magnet  or  coil 
displacement,  the  magnetic  damping  can  be  calculated. 

In  the  analysis,  various  assumptions  can  be  made  to  make  the  analysis  easier. 
The  motion  of  a  coil  or  magnet  can  be  calculated  numerically  on  the  basis  that  no 
mechanical  and  aerodynamic  damping  is  present.  From  the  time  histories  of  the 
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coil  or  magnet  motion,  magnetic  damping  can  be  calculated.  For  example,  consider 
a  single-degree-of-freedom  magnet  levitated  on  a  guideway.  From  the  oscillation 
frequency  and  time  constant  in  terms  of  the  decay  or  amplification  amplitude,  the 
magnetic  stiffness  and  damping  can  be  calculated. 

In  an  experiment,  two  or  three  series  of  tests  must  be  performed  to  obtain 
magnetic  stiffness  and  damping. 

•  Structural  stiffness  and  damping  without  magnetic  field  and 
aerodynamic  effect  can  be  measured  from  the  time  decay  of  the 
structural  response  to  an  impact. 

•  In  aerodynamic  stiffness  and  damping,  a  dummy  magnet  or  coil  can 
be  used  to  measure  the  combination  of  stiffness  and  damping  of  the 
structure  and  aerodynamic  effect.  From  the  resultant  stiffness  and 
damping,  the  structural  stiffness  and  damping  can  be  subtracted  to 
obtain  aerodynamic  stiffness  and  damping. 

•  With  a  magnetic  field,  the  total  stiffness  and  damping  will  consist  of 
magnetic,  aerodynamic,  and  structural  stiffness  and  damping. 
Magnetic  stiffness  and  damping  can  be  obtained  by  subtracting 
structural  and  aerodynamic  stiffness  and  damping. 

The  indirect  method  requires  a  series  of  tests  to  provide  the  magnetic  stiffness  and 
damping.  The  system  must  be  clearly  quantified  before  reliable  results  of  magnetic 
stiffness  and  damping  can  be  obtained. 


3.2  Direct  Method 

3.2.1  Analytical/Numerical  Method 

The  general  procedure  to  calculate  the  magnetic  forces  on  the  basis  of  unsteady 
motion  can  be  summarized  as  follows: 

Choose  a  specific  case,  such  as  a  coil  or  magnet  oscillating  over  an  aluminum 
conducting  sheet  (Fig.  1).  The  coil  is  levitated  at  height  z  and  moving  in  the  x 
direction  at  a  speed  v.  The  magnetic  force  components  acting  on  the  coil  are  fx,  fy, 
and  fz  in  the  three  directions,  x,  y,  and  z.  The  steps  described  below  can  be  taken  to 
calculate  the  needed  magnetic  forces. 
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3.2. 1 . 1  Steady  Magnetic  Force 

Set  x  =  vt,  y  =  0,  and  z  =  Zo,  calculate  static  magnetic  forces 


fxo  =  Fxo  (v,  zo,  I,  h,  a), 


fyo  =  Fyo  (V,  Zo,  I,  h,  O), 


(6) 


fzo  =  Fz0  (v,  Zo,  I,  h,  c), 


where  I  is  coil  current,  h  is  track  thickness,  and  c  is  track  conductivity. 


3.2.1 .2  Unsteady  Magnetic  Forces 

Set  the  condition  of  the  system  to  one  of  the  following  three  conditions: 

1.  Oscillations  in  the  propulsion  direction 
z  =  z0 

y  =  0  (7) 

x  =  vt  +  xi  cos  (cot) 

2.  Oscillations  in  the  guidance  direction 
z  =  z0 

y  =  y\  cos  (cot)  (8) 

x  =  vt 

3.  Oscillations  in  the  lift  direction 
Z  =  Z0  +  Z\  cos  (cot) 

y  =  0  (9) 

X  =  vt 


where  xlt  ylf  and  are  the  oscillation  amplitudes  of  the  coil  in  the  propulsion, 
guidance,  and  lift  directions,  respectively,  and  co  is  the  circular  frequency  of 
oscillations. 

Under  any  of  the  three  conditions,  the  unsteady  magnetic  force  fi  can  be 
calculated.  For  example,  when  the  coil  is  oscillating  in  the  propulsion  direction,  the 
three  components  are 
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fi  =  Fi (v,  z0,  xv  I,  to,  t,  h, a),  is  1,2, 3.  (10) 

From  Eqs.  6,  7,  and  10,  magnetic  damping  and  stiffness  matrices  can  be  calculated. 
The  magnetic  damping  and  stiffness  for  other  conditions  can  be  calculated  in  a 
similar  manner. 

It  is  believed  that  some  motion-dependent  magnetic  forces  for  simple 
geometries  can  be  predicted  analytically  or  numerically.  If  not,  they  can  be 
measured  experimentally. 


3.2.2  Experimental  Method 

An  experimental  method  has  been  developed  to  measure  the  unsteady  fluid 
forces  acting  on  a  group  of  tubes  in  cross  flow  (Chen  1987;  Chen  et  al.  10,>3;  Tanaka 
and  Takahara  1982).  The  same  technique  can  be  applied  to  measure  magnetic 
forces. 


3.2.2. 1  Quasisteady-Motion  Theory 

The  magnetic  forces  depend  on  the  deviation  from  a  reference  state,  i.e.,  the 
motion-dependent  magnetic  forces  depend  only  on  uj,  but  not  iij  and  iij,  so  that 

3 

fi  =  Xkijuj’  i,  j  =  1,2,3.  (11) 

j=l 


In  this  case,  the  magnetic  forces  are  determined  uniquely  by  coil  or  magnet  position. 
All  elements  of  magnetic  stiffness  kjj  can  be  obtained.  To  determine  kjj,  the 
magnetic-force  component  fj  is  measured  as  a  function  of  uj.  Stiffness  is  given  by 


j«L 

duj 


(12) 


3.2.2.2  Unsteady-Motion  Theory 

If  the  displacement  component  Uj  is  excited,  its  displacement  is  given  by 
uj  =  d  j  exp(V-l  tot), 


(13) 
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where  dj  is  the  excitation  amplitude  of  the  coil  or  magnet.  The  motion-dependent 
magnetic  force  of  the  component  f,  acting  on  a  vehicle  is  given  by 

fi  =  [ay  cos(i|fy)  +  ay  sin(\|/y)]d  j  exp  wt),  ( 14) 

where  ay  is  the  magnetic  force  amplitude  and  vyy  is  the  phase  angle  between  the 
magnetic  force  and  the  vehicle  displacement  Uj.  These  values  are  measured 
experimentally. 

Using  Eqs.  5  and  13,  we  can  also  write  the  motion-dependent  magnetic  force 
component  as 

fi  =  [-myco2  +  V-l  co  Cy  +  ky  jdj  exp(V-l  cot).  ( 15) 

A  comparison  of  Eqs.  14  and  15  and  solving  for  cy  and  ky  yields 
Cy  =aySin(v|/y)/cO, 

(16) 

mu=[ky-aycos(yy)]/©2. 

Based  on  Eqs.  12  and  16,  all  motion-dependent  magnetic-force  matrices  can  be 
determined  from  two  experiments:  quasisteady-motion  and  unsteady  motion. 

If  my  and  cy  are  of  no  concern,  the  experiment  that  uses  quasisteady  motion  is 
sufficient  to  determine  ky. 


4  Experimental  Setup _ 

The  general  experimental  setup,  which  includes  a  rotating  drum,  shaker,  force 
transducer,  and  displacement  transducer,  is  shown  in  Fig.  2.  The  aluminum  wheel 
is  a  rotating  drum  covered  with  a  14.61-cm  (5-3/4-in.)- wide,  6.35-mm  ( 1/4-in. )-thick 
aluminum  sheet  rolled  to  a  diameter  of  26.99  cm  (10  5/8  in.).  The  rotating  speed  can 
vary  from  0  to  3500  r/min;  therefore,  the  speed  of  the  wheel  is  0  to  50  m/s.  The 
shaker  provides  proper  excitation  force  at  given  frequencies;  and  the  impedance 
transducer  measures  the  displacement  of  the  supporting  bar. 

A  schematic  representation  of  a  force  transducer  and  magnet  support  is  shown 
in  Fig.  3.  The  magnet  is  25.4  x  50.8  x  3.18  mm.  The  magnet  is  supported  by  an 
aluminum  bar  2.54  cm  (1.0  in.)  wide,  22.86  cm  (9.0  in.)  long  (A  to  E),  and  1.27  cm 
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(0,5  in.)  thick.  The  supporting  bar  is  connected  to  the  shaker  at  one  end  and  to  the 
magnet  at  the  other.  At  one  end  it  is  machined  to  a  thickness  of  =4.78  mm 
(0.188  in.).  The  width  is  reduced  to  1.27  cm.  One  set  of  strain  gauges  is  placed  on 
the  surface  of  the  smaller  section.  The  strain  gauges  measure  the  force  components 
in  the  lift  direction  with  a  sensitivity  of  =1  v  for  120  g  of  force  acting  on  the  middle 
of  the  magnet. 

The  supporting  bar  with  the  magnet  fixed  to  the  shaker  at  the  upper  end  is 
fairly  rigid  with  a  natural  frequency  of  >100  Hz.  In  the  tests,  the  frequencies  of 
interest  are  <20  Hz.  Therefore,  the  response  of  the  supporting  bar  will  not  affect 
the  measurement  of  magnetic  forces. 

The  force  transducer  is  calibrated  by  both  static  and  dynamic  methods.  In  the 
static  method,  the  supporting  bar  is  attached  to  the  shaker  at  one  end  and  a  given 
force  is  applied  to  the  middle  of  the  active  length  (DE).  Because  the  supporting  bar 
is  an  elastic  beam,  the  applied  force  can  be  correlated  with  the  displacement 
measured  by  the  displacement  transducer  and  the  change  of  the  clearance  between 
the  drum  and  magnet.  Therefore,  once  the  displacement  is  measured  by  the 
displacement  transducer,  the  steady  magnetic  force  and  the  clearance  between  the 
magnet  and  drum  can  be  calculated.  The  steady  magnetic  lift  force  f  (g)  is  given  by 

f=  403.9  d,  (17) 

where  d  (mm)  is  the  displacement  measured  by  the  displacement  transducer.  The 
change  of  the  clearance  G  (mm)  between  the  magnet  and  drum  is  given  by 

G  =  0.004117  f.  (18) 

This  means  that  for  a  static  force  applied  to  the  magnet,  the  clearance  will  be 
increased  by  G. 

In  the  dynamic  method,  the  force  transducer  is  excited  at  a  given  frequency 
and  amplitude.  Then,  the  inertia  forces  due  to  the  sinusoidal  oscillations  are  used 
to  determine  the  calibration  constant. 

All  calibration  constants  used  in  the  measurement  of  motion-dependent 
magnetic  forces  are  determined  by  the  dynamic  method.  The  inertia  force  is  given 
by 


fj  =  mo^d  cos(o)t), 


(19) 
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where  G)  is  the  oscillation  frequency  (rad/sec),  d  is  oscillation  peak  amplitude,  and  m 
is  the  mass  of  the  active  element,  which  is  shown  in  Fig.  3b.  The  active  element 
includes  the  lower  portion  of  the  supporting  bar  (from  the  middle  of  the  strain  gauge 
to  the  lower  end),  magnet,  and  aluminum  or  copper  plate.  In  testing,  the  inertia 
force  fj  is  measured  for  a  given  excitation  0)  and  vibration  amplitude  d.  In  the 
calibration,  to  varies  from  =2  to  8  Hz  and  d  varies  from  0.8  to  1.1  mm.  From  Eq.  19, 
the  ratio  of  fj/d  is  proportional  to  to2.  Figure  4  shows  fj/d  as  a  function  of  excitation 
frequency  for  a  typical  case,  Test  A.  Theoretically,  the  ratio  fj/d  should  be 
proportional  to  to2;  the  actual  power  is  very  close  to  2.0. 


5  Test  Procedures  and  Data  Analysis _ 

A  flow  diagram  of  the  instrumentation  and  shaker  (electrodynamic  exciter)  is 
shown  in  Fig.  5.  The  shaker  provides  a  sinusoidal  displacement  with  a  frequency 
varying  from  -2  to  15  Hz.  Displacement  and  force  signals  can  be  filtered  first  to 
eliminate  low-  and  high-frequency  noises.  These  signals  are  analyzed  to  obtain  the 
oscillation  displacement  of  the  supporting  bar,  the  magnitudes  of  the  forces  acting 
on  the  magnet,  and  the  phase  between  the  motion-dependent  magnetic  force  and 
the  displacement  of  the  supporting  bar. 

Three  series  of  tests  were  performed,  as  shown  in  Fig.  3b  and  Table  1. 

A.  Magnet:  Only  the  magnet  is  attached  to  the  force  transducer.  The 
clearance  between  the  magnet  and  rotating  drum  can  be  set  at 
different  values. 

B.  Magnet  and  Aluminum  Plate:  An  aluminum  plate  4.75  mm  thick, 

7.62  cm  long,  and  5.08  cm  wide  is  attached  to  the  magnet.  The 
clearance  between  the  magnet  and  rotating  drum  was  set  at 
10.5  mm. 

C.  Magnet  and  Copper  Plate:  A  copper  plate  1.55  mm  thick,  10.16  cm 
long,  and  7.62  cm  wide  is  attached  to  the  magnet.  The  clearance 
between  the  magnet  and  rotating  drum  was  set  at  10.5  mm,  while  the 
copper  plate  can  be  placed  in  the  gap  with  different  clearance  with 
respect  to  the  rotating  drum. 

The  magnet  may  be  replaced  by  a  dummy  steel  plate  with  the  same  s'ze  and 
same  weight  as  the  magnet.  With  the  dummy  plate,  the  dynamic  characteristics  of 
the  supporting  structure,  as  well  as  the  aerodynamics,  can  be  characterized.  In  all 
of  the  experiments,  two  tests  must  be  performed:  one  with  the  magnet  and  the  other 
with  the  dummy  steel  plate. 


11 


Table  1.  Summary  of  test  cases 


Test  Cases 

Material 

Clearance,  mma 

Total  Mass,  gb 

A.l 

Magnet  Only 

10.5 

96.2 

A.2 

Magnet  Only 

7.5 

96.2 

A.3 

Magnet  Only 

6.5,  8.5,  10.5, 

12.5,  and  14.5 

96.2 

B 

Magnet  and 
Aluminum  Plate 

4.5 

133 

C.l 

Magnet  and 
Copper  Plate 

7.7 

186.5 

C.2 

Magnet  and 
Copper  Plate 

3 

226 

aGap  between  the  aluminum  sheet  of  the  drum  and  the  magnet,  aluminum 
plate,  or  copper  plate. 

bMass  of  active  element  shown  in  Fig.  3b,  including  supporting  aluminum 
bar  below  middle  of  the  strain  gauge,  magnet,  and  aluminum  or  copper 
plate. 


In  the  three  series  of  tests,  the  clearance  G  may  vary.  The  test  cases  are 
summarized  in  Table  1.  Test  Case  A  is  the  reference  case,  whereas  the  others  are 
performed  to  enhance  our  understanding  of  the  effects  of  the  additional  plate  on  the 
motion-dependent  magnetic  forces. 


6  Test  Results 


The  steady  magnetic  force  depends  on  the  rotating  speed  of  the  drum  and  the 
clearance.  Tests  were  performed  for  clearances  from  =6.5  to  14.5  mm  for  the 
different  cases. 

The  measured  dynamic  force  is  given  as  follows: 
dt2  dt 


(20) 
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where  u  is  the  displacement,  m  is  the  mass  of  the  active  element,  denoted  by  total 
mass  in  Table  1,  and  C  and  K  are  damping  coefficient  and  stiffness  to  be 
determined,  respectively.  Note  that  m  should  also  include  magnetic  mass,  as  shown 
in  Eq.  5.  However,  because  the  oscillation  frequency  of  maglev  is  fairly  low,  the 
magnetic  mass  will  be  very  small  and  is  negligible  (Iwamoto  et  al.  1974). 

From  the  time  signals  of  the  force  transducer  and  displacement,  the  RMS 
values  of  the  displacement  and  force  are  d0  and  f0,  respectively,  and  the  phase  angle 
between  the  two  is  <j>.  C  and  K  are  given  by 

C  =  f0  sin  (<}>)  /cflu0 

(21) 

K  =  mco2  -  (fo/u0)  cos  (4>). 

For  a  given  magnet  or  coil  on  a  sheet  guideway,  the  magnitude  of  the  force 
depends  on  excitation  amplitude,  frequency,  and  speed.  Figure  6a  shows  the 
magnetic  force  and  displacement  at  2  and  14  Hz  of  excitation  at  4.8  m/s,  and  Fig.  6b 
shows  magnetic  force  with  an  excitation  frequency  of  6  Hz  for  two  speeds,  0  and 
16.55  m/s.  The  signals  in  Fig.  6  were  not  filtered.  The  frequency  of  magnetic  force 
is  the  same  as  the  excitation  frequency;  however,  the  magnetic  force  is  not  in  phase 
with  the  displacement.  Depending  on  the  phase  between  the  magnet  displacement 
and  magnetic  forces,  the  magnetic  forces  may  be  an  excitation  mechanism  or  energy 
dissipation  mechanism.  Figures  6a  and  b  also  show  the  loci  of  the  unsteady 
magnetic  force  with  displacement.  In  Fig.  6a,  at  2  Hz,  and  in  Fig.  6b,  at  16.55  m/s, 
the  magnetic  force  is  negative  when  the  displacement  is  positive  and  vice  versa; 
therefore,  it  acts  as  an  excitation  force.  In  Fig.  6a,  at  14  Hz,  and  in  Fig.  6b,  at  0  m/s, 
the  magnetic  force  is  positive  when  the  displacement  is  positive;  therefore,  it  acts  as 
a  damping  force.  A  damping  force  will  contribute  to  energy  dissipation,  which  will 
reduce  maglev  response,  whereas  an  excitation  force  will  contribute  to  negative 
damping,  which  may  result  in  dynamic  instability. 

The  phase  angle  is  defined  as  the  angle  between  the  magnetic  force  and  the 
magnetic  displacement.  This  angle  can  be  obtained  from  the  correlation  between 
the  two  time  histories  of  magnetic  force  and  displacement,  or  from  the  time  history 
plots  given  in  Fig.  6. 


6.1  Test  A.1 

Figures  7  and  8  show  the  steady  magnetic  force,  RMS  displacement,  RMS 
magnetic  force,  phase  angle,  and  calculated  damping  coefficient  and  stiffness  as  a 
function  of  rotating  speed  at  two  excitation  frequencies,  2  Hz  and  6  Hz.  Results  for 
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both  magnet  and  dummy  are  given.  The  net  magnet  stiffness  and  damping  are 
given  in  Fig.  9.  Once  the  drum  is  rotating,  a  steady  lift  force  is  imposed  on  the 
magnet  and  the  clearance  increases.  The  change  of  the  clearance  can  be  calculated 
from  Eq.  18.  The  structural  damping  does  not  change  with  the  clearance,  whereas 
the  aerodynamic  force  is  small,  when  compared  with  magnetic  force,  and  the  effect 
of  the  clearance  due  to  static  deflection  of  the  force  transducer  is  also  small. 
Therefore,  when  the  dummy  was  installed,  the  clearance  was  not  adjusted  to 
account  for  the  effect  of  the  small  change  of  clearance  due  to  static  deflection  on  the 
aerodynamic  force. 

The  following  general  characteristics  were  noted: 

•  The  steady  magnetic  force  increases  with  speed.  The  results  compare 
reasonably  well  with  previous  measurements  (Cai  et  al.  1992). 

•  The  phase  angle  between  the  force  and  displacement  change 
drastically  at  the  characteristic  speed  (about  7  m/s),  which  is  equal  to 
2/pah  [p  is  the  permeability,  a  is  the  aluminum  sheet  conductivity, 
and  h  is  the  smaller  of  aluminum  sheet  thickness  and  skin  thickness 
(Sinha  1987)];  therefore  the  damping  coefficient  changes  sign  when 
the  speed  is  increased  to  the  characteristic  speed. 

•  Damping  and  stiffness  do  not  change  significantly  with  excitation 
frequency. 

•  The  aerodynamic  damping  and  stiffness  are  fairly  small. 

Figure  10  shows  damping  coefficient  and  stiffness  as  a  function  of  excitation 
frequency  for  three  speeds.  Also  noted  in  Figs.  7-9,  the  stiffness  and  damping  are 
relatively  independent  of  the  excitation  frequency. 

The  effect  of  the  excitation  amplitude  on  the  magnetic  damping  and  stiffness 
was  also  tested.  The  results  are  given  in  Fig.  11  at  33.6  m/s  and  four  frequencies,  2, 
6,  10,  and  14  Hz.  In  the  amplitude  range  tested,  the  effect  is  very  small.  Therefore, 
for  small  oscillations,  the  magnetic  damping  and  stiffness  can  be  considered  as 
independent  of  oscillation  amplitude. 


6.2  TestA.2 

Figures  12  and  13  show  the  results  for  2  Hz  and  6  Hz.  The  aerodynamic  and 
structural  contributions  are  subtracted,  based  on  the  results  of  10.5-mm  clearance, 
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because  they  are  insensitive  to  spacing  in  this  range  of  clearance.  Magnetic 
damping  and  stiffness  as  a  function  of  speed  are  given  in  Fig.  14. 

The  difference  between  Test  A.1  and  Test  A.2  is  the  clearance.  In  Test  A.l,  the 
clearance  at  0  speed  is  10.5  mm,  whereas  in  Test  A.2,  it  is  6.5  mm.  The  general 
trend  of  magnetic  damping  and  stiffness  are  similar.  At  small  clearance,  the 
magnitude  of  magnetic  damping  and  stiffness  are  larger  at  smaller  clearances.  It  is 
also  noted  that  at  zero  speed,  magnetic  damping  is  larger  with  a  smaller  clearance. 
This  agrees  with  the  experimental  data  of  Moon  (1977). 


6.3  Test  A.3 

To  show  the  effect  of  the  clearance  and  excitation  frequency,  tests  were 
performed  at  50  m/s  for  clearances  of  6.5,  8.5,  10.5,  12.5,  and  14.5  mm  and 
frequencies  of  2-6  Hz.  Figures  15  and  16  show  the  data  obtained  when  testing  with 
magnet  and  dummy.  In  these  two  figures,  the  displacement  is  referred  to  as  static 
displacement.  With  the  dummy,  the  displacement  is  zero  because  there  is  no 
magnetic  force.  Figure  16  shows  that  the  damping  and  stiffness  are  practically  zero 
at  various  frequencies. 

Figures  17  and  18  show  damping  and  stiffness  as  a  function  of  excitation 
frequencies  for  five  clearances.  Both  the  damping  and  stiffness  coefficients  are 
relatively  independent  of  the  excitation  frequency.  The  absolute  magnitudes  of 
magnetic  stiffness  and  damping  increase  with  reducing  clearance.  At  50  m/s, 
magnetic  damping  is  negative. 


6.4  Test  B 

In  this  test,  an  aluminum  plate  (76.2  x  50.8  x  4.75  mm)  is  attached  to  the 
magnet  with  a  clearance  of  1.3  mm  between  the  surface  of  the  aluminum  plate  and 
the  magnet  (Fig.  3b). 

The  results  for  magnetic  damping  and  stiffness  are  given  in  Figs.  19-21.  To 
verify  the  repeatability,  two  tests  were  performed  to  measure  the  data  as  a  function 
of  excitation  frequency.  The  results  are  given  in  Fig.  22.  The  difference  is  very 
small.  This  demonstrates  the  repeatability  of  the  test  setup. 

The  effect  of  the  additional  aluminum  on  the  damping  and  stiffness  can  be 
examined  by  comparing  the  results  with  Test  A.l.  When  comparing  Figs.  7a  and 
19a,  Figs.  8a  and  20a,  Figs.  9  and  21,  we  observe  that  the  difference  is  very  small. 
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This  means  that  in  this  case  the  additional  plate  attached  to  the  magnet  does  not 
change  the  magnetic  damping  and  magnetic  stiffness. 


6.5  Test  C.1 

A  copper  plate  is  attached  to  the  magnet.  Figures  23-25  show  the  results  for 
case  C.l.  The  effect  of  the  additional  copper  plate  on  the  damping  and  stiffness  can 
be  examined  by  comparing  the  results  with  Test  A.l.  When  comparing  Figs.  7a  and 
23a,  Figs.  8a  and  24a,  Figs.  9  and  25,  we  see  that  the  difference  is  very  small.  This 
means  that  the  additional  copper  plate  attached  to  the  magnet  does  not  change  the 
magnetic  damping  and  magnetic  stiffness  significantly.  Figures  20a  and  24a  show 
that  the  dynamic  magnetic  force  at  zero  speed  is  larger  than  at  5  m/s. 


6.6  TestC.2 

In  this  test,  the  copper  plate  is  moved  closer  to  the  drum  while  the  magnet  is  at 
the  same  location.  Figures  26-28  show  Test  C.2  data  measured  at  2  Hz  and  6  Hz. 
The  results  for  magnetic  damping  and  stiffness  are  given  in  Fig.  28. 

Figures  29a-e  show  the  data  for  the  magnet  with  the  copper  plate  as  a  function 
of  excitation  frequency  at  five  speeds.  The  corresponding  data  for  the  dummy  with 
copper  plate  are  given  in  Fig.  30.  The  magnetic  damping  and  stiffness  are  given  in 
Fig.  31. 

A  comparison  of  Figs.  25  and  28  shows  that  the  magnetic  damping  increases 
slightly  once  the  speed  is  not  zero.  This  means  that  when  the  copper  plate  is  closer 
to  the  drum,  the  magnetic  damping  will  increase. 

Figure  31  shows  magnetic  damping  and  stiffness  as  a  function  of  excitation 
frequency  for  a  series  of  speeds.  At  0  and  5  m/s  the  magnetic  damping  is  positive. 
At  three  other  speeds,  the  magnetic  damping  values  are  approximately  the  same 
and  decrease  slightly  with  frequency.  The  magnetic  stiffness  is  zero  at  0  m/s  and 
varies  from  =370  N/m  to  470  N/m  at  5  m/s  as  the  frequency  is  varied  from  2  to  6  Hz. 

On  the  basis  of  the  experimental  data,  some  general  remarks  can  be  made. 

•  Damping,  including  support  structures  and  aerodynamics,  is  fairly 
small  in  the  test  setup.  At  various  speeds,  the  damping  values 
without  a  magnet  are  positive,  in  general.  In  some  cases,  because  the 
damping  values  are  small,  some  error  in  the  measurements  results  in 
very  small  negative  values.  The  small  errors  in  the  direct  method  of 


measurement  do  not  affect  the  general  characteristics  of  the  magnetic 
forces. 

Without  a  magnet,  the  stiffness  is  practically  zero.  This  implies  that 
the  stiffness  due  to  aerodynamic  forces  is  practically  zero. 

At  zero  speed,  the  magnetic  damping  is  always  positive.  As  the 
clearance  between  the  magnet  and  drum  is  reduced,  the  damping 
coefficient  increases.  The  stiffness  is  zero  at  zero  speed. 

When  the  speed  is  lower  than  the  characteristic  speed,  magnetic 
damping  is  positive.  As  the  speed  increases,  it  becomes  negative. 

At  high  speed,  which  is  much  greater  than  the  characteristic  speed, 
the  damping  and  stiffness  are  almost  independent  of  the  speed;  but 
they  depend  on  the  clearance  between  the  magnet  and  drum.  As  the 
clearance  is  reduced,  the  magnitude  of  magnetic  damping  and 
stiffness  increase.  The  magnetic  damping  is  always  negative  with 
the  absolute  magnitude  increasing  with  decreasing  clearance. 

In  the  parameter  range  tested,  the  effect  of  excitation  frequency  on 
the  magnetic  damping  and  stiffness  is  very  small.  This  implies  that 
the  data  measured  at  a  specific  frequency  can  be  used  for  other 
frequencies  without  much  error. 

The  effect  of  excitation  amplitude  on  the  magnetic  damping  and 
stiffness  is  also  very  small.  However,  if  the  motion  is  very  large, 
nonlinear  effects  are  expected  to  become  important.  If  one  is 
interested  in  the  initiation  of  large  motion,  the  linear  theory  is 
applicable  and  the  magnetic  forces  measured  for  small  amplitudes 
will  be  applicable.  If  one  is  interested  in  the  dynamics  response  with 
large-amplitude  oscillations,  such  as  oscillations  in  the  instability 
region,  nonlinear  effects  must  be  included. 

The  passive  damping  that  results  from  the  aluminum  and  copper 
plates  in  the  test  is  small.  Other  geometrical  shapes  or  different 
materials  at  various  temperatures  will  be  needed  to  generate 
sufficient  passive  damping  for  maglev. 

With  a  fixed  clearance  between  the  magnet  and  rotating  drum,  the 
magnetic  damping  due  to  the  copper  plate  placed  between  the 
magnet  and  drum  increases  with  the  decrease  of  the  clearance 
between  the  plate  and  drum. 
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The  general  magnetic  damping  characteristics  are  similar  to  the  train  and 
track  loop  (Iwamoto  et  al.  1974).  Magnetic  damping  is  negative  once  the  speed  is 
larger  than  the  characteristic  speed. 


7  Applications _ 

In  this  study,  we  measured  the  lift  force  due  to  the  oscillation  in  the  lift 
direction.  The  magnetic  damping  and  stiffness  are  due  to  the  motion  of  itself;  in 
general,  they  are  called  self-induced  damping  and  stiffness.  In  reality,  once  the 
magnet  is  given  a  motion  in  one  direction,  it  will  also  induce  magnetic  forces  in 
other  directions.  Specifically,  if  the  magnet  is  excited  in  the  lift  direction,  three 
components  of  the  induced  forces  in  the  lift,  drag,  and  guidance  directions  will  be 
induced.  In  experiments,  similar  designs  of  the  support  and  force  transducers  can 
be  developed  to  measure  the  three  force  components  that  are  attributable  to  an 
excitation  in  a  direction.  From  the  time  histories  of  these  three  components  of 
magnetic  forces,  as  well  as  the  phases  between  the  force  components  and 
displacement,  the  magnetic  damping  and  stiffness  in  the  three  directions  can  be 
calculated.  The  magnetic  damping  and  stiffness  in  one  direction  due  to  motion  in 
another  direction  are  called  coupled  magnetic  damping  and  stiffness.  For  example, 
for  a  magnet  with  the  force  in  three  directions,  the  resulting  magnetic  damping  and 
stiffness  are  the  same  as  those  given  in  Eq.  5. 

The  method  presented  in  this  report  is  applicable  to  general  cases.  If  the 
magnetic  damping  and  stiffness  cannot  be  calculated,  the  direct  method  presented 
here  will  be  useful.  In  some  cases  in  which  the  coupled  damping  and  stiffness  are 
not  important,  the  indirect  method  may  be  used  to  measure  the  self-induced 
damping  and  stiffness.  When  the  coupling  effect  is  important,  it  will  be  difficult  to 
identify  all  of  the  elements  of  magnetic  damping  and  stiffness  matrices  from  the 
results  of  the  indirect  method. 

The  experimental  data  presented  in  Section  6  can  be  applied  to  the  dynamics  of 
a  maglev  vehicle.  Once  the  magnetic  damping  and  stiffness  are  known,  the 
response  and  stability  of  maglev  systems  can  be  predicted.  Because  we  measured 
only  the  lift  force  in  this  study,  only  the  oscillation  in  the  lift  direction  is  considered. 

Consider  a  maglev  vehicle,  with  the  mass  m,  rising  in  a  guideway.  The 
applicable  equation  is 

m  +  (Cs  +C)^  +  Ku  =  q(t), 

dt2  dt 


(22) 
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where  m  is  the  total  mass  of  the  vehicle,  Cs  is  structural  damping,  including 
aerodynamic  damping;  C  is  magnetic  damping;  K  is  magnetic  stiffness;  and  q(t)  is 
external  excitation. 


The  natural  frequency  and  modal  damping  ratio  are 


f  = 


2k 


C s+C 

^  2VKm' 


(23) 


For  a  given  vehicle,  C  and  K  depend  on  clearance  and  speed;  both  natural  frequency 
and  modal  damping  will  depend  on  clearance  and  speed. 

When  the  damping  value  is  zero,  i.e. 


C  =  o 


or 


(24) 


Cs  +  C  =  0, 


instability  will  occur. 

Consider  a  specific  case.  A  vehicle  with  mass  equal  to  150  g  travels  at  a  speed 
of  50  m/s.  The  magnetic  damping  and  stiffness  are  given  in  Fig.  17.  The  natural 
frequency  and  modal  damping  due  to  magnetic  forces  can  be  calculated  from  Eq.  23 
assuming  Cs  =  0.  The  results  are  given  in  Fig.  32.  As  the  clearance  increases,  the 
natural  frequency  and  negative  magnetic  damping  decrease. 

To  predict  the  dynamic  response  of  a  maglev  system,  magnetic  damping  and 
stiffness  matrices  must  be  quantified.  Using  the  direct  method,  either  analytically 
or  experimentally,  to  characterize  magnetic  damping  and  stiffness,  including  self- 
induced  and  coupling  terms,  we  can  predict  a  detailed  analysis  of  maglev  stability 
and  response. 


8  Closing  Remarks 


Motion-dependent  magnetic  forces  are  the  controlling  elements  in  the 
prediction  of  dynamic  response  and  stability  of  maglev  systems.  At  this  time,  very 
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limited  studies  have  been  published  and  many  stability  issues  remain  unresolved. 
An  integrated  analytical/experimental  study  has  been  initiated  at  Argonne  National 
Laboratory  (ANL)  to  develop  an  experimental  method  and  computational  tools  for 
the  study  of  motion-dependent  magnetic  forces.  This  report  presents  a  general 
analytical  approach  and  an  experimental  technique  to  measure  the  motion- 
dependent  magnetic  forces  by  a  direct  method.  The  method  is  applicable  to  general 
cases  and  will  provide  the  necessary  information  on  the  magnetic  forces  of  rnaglev 
systems  and  resolve  the  issues  of  maglcv  stability. 

Application  of  the  experimental  method  presented  in  this  report  can 

•  Provide  the  key  elements  of  magnetic  forces  for  stability  analysis  of 
maglev  systems.  Once  the  magnetic  damping  and  stiffness 
coefficients  for  the  coil  and  magnet  are  known,  these  coefficients  can 
be  applied  to  maglev  systems. 

•  Characterize  the  general  characteristics  of  magnetic  damping  and 
magnetic  stiffness.  It  might  be  possible  to  identify  some  nondimen- 
sional  parameters  that  can  be  used  to  determine  the  role  of 
quasisteady-  and  unsteady-motion  theories. 

•  Provide  the  answer  to  the  question:  what  are  the  dominant  elements 
that  cause  instability  in  maglev  systems?  Magnetic  damping  or 
magnetic  stiffness?  Self-induced  magnetic  damping  or  coupling  of 
magnetic  stiffness? 

•  Enable  the  development  of  a  more  efficient  method  of  calculating 
magnetic  damping  and  stiffness  coefficients  by  enhancing  our 
understanding  of  the  general  characteristics  of  magnetic  damping 
and  stiffness. 

•  Enable  us  to  perform  a  complete  stability  analysis  of  maglev  systems 
by  exploring  motion-dependent  magnetic  forces.  It  will  be  possible  to 
resolve  a  series  of  unsettled  issues  on  maglev  stability. 

This  study  represents  the  first  step  in  the  development  of  magnetic  damping 
and  stiffness  for  maglev.  The  experimental  technique  that  uses  the  direct  method 
will  provide  the  necessary  elements  of  magnetic  forces.  As  a  follow-up  investigation, 
two  series  of  tests  will  be  performed  at  ANL:  (1)  conducting  plates  or  coils  will  be 
inserted  between  the  magnet  and  drum  to  study  the  role  of  passive  damping;  and 
(2)  a  magnet  will  be  excited  in  one  direction  and  the  forces  in  several  directions  will 
be  measured  simultaneously.  The  results  from  these  studies  will  improve  our 
understanding  of  the  dynamics  of  maglev  systems. 
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Fig.  2.  Experimental  setup;  (a)  schematic 
representation  and  (b)  photograph 
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Fig.  3.  Force  transducer  and  magnet 
support;  (a)  overall  view;  (b)  active 
elements  for  three  test  cases 
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Fig. 


Force-displacement  ratio  as  a  function  of  excitation  frequency  for  Test  A 


26 


Fig.  5.  Flow  diagram  of  data  analysis  instrumentation 
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Fig.  11.  (a)  Magnetic  damping  and  (b)  magnetic  stiffness  as  a  function  of  excitation 
amplitude,  Test  A.  1 
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Fig.  14.  Magnetic  damping  and  stiffness.  Test  A.2 
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Fig.  18.  Magnetic  stiffness  as  a  function  of  frequency  for  various 
clearances.  Test  A.  3 
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Fig.  32.  Natural  frequency  and  magnetic  modal  damping  ratio 
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